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Abstract 

Sensitive imaging X-ray observations of young stellar clusters (YSCs, 
ages < 10 Myr) are valuable tools for the acquisition of an unbiased cen- 
sus of cluster members needed for Initial Mass Function (IMF) studies. 
Several dozen YSCs, both nearby and across the Galactic disk, have 
been observed with the Chandra and XMM-Newton satellites, detect- 
ing > 10, 000 low-mass cluster members. Many of these samples should 
be nearly complete down to ~ 1 Mq. An important additional benefit 
is that the YSC X-ray luminosity function appears to be universal with 
a lognormal shape, providing a new standard candle for measurement 
of YSC distances. 



1. Measuring IMFs in young stellar clusters 

Measuring the stellar Initial Mass Function (IMF) requires complete 
and unbiased samples which are often sought in young stellar clusters 
(YSCs) that have not undergone significant dynamical evolution. We 
consider here YSCs with ages < 10 Myr where most of the lower mass 
stars are on the convective Hayashi pre-main sequence (PMS) evolution- 
ary tracks. 

A major difficulty in obtaining a reliable census of YSCs is that the few 
lying at distances d < 400 pc are spatially extended and not strongly con- 
centrated. The richest is the Sco-Cen OB Association with several thou- 
sand members, but it subtends several thousand square degrees in the 
southern sky. Others like the Corona Australis, Perseus or Chamaeleon 
clouds have dozens of members and typically subtend > 1 square de- 
gree. Individual Taurus-Auriga clouds or Bok globules may be smaller, 
but their stellar populations are too poor for IMF studies. 
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These large angular sizes require that cluster members of nearby 
clouds must be efficiently discriminated from the large, often overwhelm- 
ing, number of unrelated foreground and background stars. This dis- 
crimination is traditionally achieved using surveys for stars with Ha 
emission and/or near-infrared photometric excesses. These methods ef- 
ficiently select protostars and 'classical T Tauri' (CTT) stars, but often 
undersample the large 'weak-lined T Tauri' (WTT) population (Feigel- 
son & Montmerle 1999). As the mass-dependency of disk evolution is 
not well known, correction to the full PMS population for IMF study 
is uncertain. Analysis of the K-hand star counts towards YSCs can 
give a more complete census but in an indirect fashion (Lada &: Lada 
1995). It often requires major corrections for unrelated stars and, as 
individual members are not individually identified, the IMF is derived 
by a model-dependent conversion of the i^-band distributions to mass 
distributions. 

IMF studies may be advantageous in more distant YSCs, but here 
additional problems present themselves. The members are fainter than 
in closer clusters, so that large telescopes with low-noise detectors be- 
come necessary to detect the lowest mass cluster members. As distance 
increases, it becomes increasingly difficult to resolve multiple systems. 
Both observational and theoretical studies suggest that most stars form 
in binary or multiple systems. A complete census is also hindered by 
the wide range of extinctions often exhibited by YSC members. 

In light of these constraints, one YSC has unique merits for IMF 
studies: the Orion Nebula Cluster (ONC). The ONC lies at d ~ 450 pc, 
sufficiently close that the best available near-infrared instrumentation 
can detect objects down to a few Jupiter masses. With ~ 2000 members 
concentrated in ~ 0.1 square degrees, it is easily studied and rich enough 
to populate the IMF up to ~ 45 Mq. Unlike other YSCs in the Orion 
molecular cloud complex, it lies on the near side of the cloud and has 
evacuated most of the intervening molecular material, so line-of-sight 
absorptions are low. While the ONC census suffers some difficulties — 
mild contamination from unrelated young stars in the background cloud 
and from foreground dispersed young clusters, incomplete enumeration 
of multiple systems — it provides a standard for stellar IMF studies that 
is far above that achieved for other YSCs. 

2. The role of X-ray surveys 

X-ray emission from PMS was initially predicted from shocks associ- 
ated with their stellar wind, but studies with the Einstein and ROSAT 
observatories indicated an origin more closely associated with solar-like 
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magnetic activity (Feigelson &: Montmerle 1999). The emission is 1 — 4 
orders of magnitude stronger than typical main sequence levels, ex- 
hibits high-amplitude flares and hotter plasmas than expected from wind 
shocks, and is mostly uncorrelated with the presence or absence of an 
infrared-excess disk. X-ray surveys discovered many WTTs whose pop- 
ulation is comparable to the CTT population even in YSCs associated 
with active star forming clouds. Although useful, surveys with these 
early satellites were limited by detector technology: low-resolution gas 
proportional counters or low-efficiency solid state microchannel plates. 

The Chandra X-ray Observatory and XMM-Newton missions are much 
better adapted to YSC studies with their high-efficiency low-noise CCD 
detectors. Chandra is particularly useful with its high-precision mirrors 
giving arcsccond imaging capability, though its field of view is limited 
to 0.08 square degrees. Together, these telescopes have imaged several 
dozen YSC populations across the Galactic disk. These include the p 
Ophiuchi, Chamaeleon, CrA, Mon R2, L1551, L1448, Serpens, Cep B 
and Sgr B2 clouds; populations associated with the Orion, Trifid, Rosette 
and Carina HII regions; the Cyg 0B2 association, Cha cluster, and 
various isolated Herbig Ae/Be stars with their companions; star forming 
regions NGC 1333, IC 348, M 8, M 16, M 17, NGC 281, NGC 1579, 
NGC 2078, NGC 2024, NGC 2264, IC 1396, RCW 38, RCW 49, RCW 
108, NGC 6334, NGC 6530, NGC 6383, NGC 3603, W 3, W 1, W 49, W 
51, and IRAS 19410-^2336; several Galactic Center YSCs; and 30 Dor in 
the Large Magellanic Cloud. The most comprehensive study underway 
is the Chandra Orion Ultradeep Project (COUP) based on a nearly- 
continuous 10-day pointing towards the ONC and embedded sources in 
OMC-1. An introduction to COUP, and references for the other regions 
published through mid-2003, are given in Feigelson (2003a). Figure 1 
shows a portion of the COUP image, and some preliminary COUP re- 
sults are presented here. 

X-ray surveys are subject to selection effects which must be carefully 
considered in the effort to achieve a well-defined YSC census. Most im- 
portantly, PMS X-ray luminosities are strongly correlated with a tangle 
of interrelated stellar properties: bolometric luminosity, mass, stellar 
surface area and volume. The — Lhoi correlation has been known for 
many years, but the linkage with other variables only emerged clearly 
in pre-COUP studies of the ONC (Flaccomio et al. 2003a; Feigelson et 
al. 2003b). Though quantitative analysis of the scatter in these rela- 
tions by the COUP has not yet been completed, we can roughly say 
that a Chandra observation of a nearby YSC with limiting sensitivity 
logL^(lim) ~ 28.0 erg/s (0.5 - 8 keV band) will detect > 90% of PMS 
stars with masses M > 0.2 Mq, and an observation of a more distant 
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Figure 1. The central 4' x 4' region of the 17' x 17' Chandra ACTS image of the Orion 
Nebula Cluster, obtained in a ~ 10-day exposure in early 2003 (Getman et al. 2004). 
Moderate-resolution spectra over the 0.5—8 keV band, giving an independent measure 
of line-of-sight absorption, and temporal information on variability are available for 
each of the 1616 COUP sources. This datasct is the foundation of a wide range of 
studies comprising the Chandra Orion Ultradeep Project (COUP). 

YSC with logLa;(lim) ~ 29.5 erg/s will detect > 90% of PMS stars with 
M > 0.8 Mq (Preibisch et al., in preparation).^ 

This result immediately indicates the power of Chandra studies for im- 
proving our knowledge of stellar populations in a wide variety of YSCs, 
and thus studying the IMF in a wide variety of conditions. In most 
YSCs beyond d ~ 1 kpc, cluster membership is currently limited to a 
small number of bright spectroscopically-confirmed OB stars and lower 
mass stars with K-band excesses. Consider, for example, the cluster illu- 
minating the HII region M 17 and its the surrounding molecular cloud. 



second second selection effect in PMS X-ray luminosities must be considered: the sub- 
population of CTT stars tends to be a factor of 2 — 3 weaker the subpopulation of WTT stars 
leading to, opposite to traditional methods, a selection bias against accreting stars (Flac- 
comio et al. 2003b). But this effect is relatively small compared to the 4 orders of magnitude 
range in PMS X-ray luminosities. 
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Only a few dozen have optical spectra and while >20,000 have JHK 
measurements, most of these stars are unrelated to the cloud (Hanson 
et al. 1997; Jiang et al. 2002). A Chandra image with logL^(lim) = 29.7 
erg/s shows 877 sources nearly all of which are cluster members (Getman 
et al., in preparation). Such X-ray selected samples should be > 90% 
complete above a well-defined mass limit, usually around 0.5 — 1.5 Mq 
for the more distant YSCs under study. 

Analysis of these fields is difficult with hundreds of faint X-ray stars 
often embedded in structured diffuse emission from large-scale OB wind 
shocks (Townsley ct al. 2003). But sophisticated data analysis methods 
have been developed (Getman et al. 2004), careful studies are underway 
and within a few years > 10,000 Chandra-discovered YSCs members 
should be published. Followup optical/near- infrared photometry and 
spectroscopy will be needed to place these stars on the HR diagram to 
estimate masses for IMF analysis. Chandra should thus provide an im- 
portant boost to comparative IMF studies for many of the YSC clusters 
listed above. 

We caution that it is still unclear whether even the most sensitive 
X-ray surveys can efficiently detect substellar PMS objects which will 
evolve into L- and T-typc brown dwarfs. Two YSC studies with sensitiv- 
ities of logL^(lim) ~ 27.0 erg/s performed to date give inconsistent re- 
sults. A deep Chandra image of the Chamaeleon I North cloud found all 
27 known cloud members including 3 probable substellar objects (Feigel- 
son & Lawson 2004). This suggests the X-ray census is complete and 
gives an IMF deficient in stars with M < 0.3 Mq compared to the ONC 
or field star IMFs. But this may be due to mass segregation favoring 
higher mass stars in the small 0.6 pc^ region covered by the Chandra 
imager, rather than intrinsically different IMF. The second study with 
logLa;(lim) ~ 27.0 erg/s is the COUP observation of the ONC. Here, 
most of the spectroscopically confirmed brown dwarfs (Slesnick et al. 
2004) are not detected, indicating that X-ray surveys will probably be 
incomplete in this low mass regime. 

3. A new distance estimator for YSCs 

Techniques for measuring distances to YSCs and their natal molecular 

clouds have hardly changed in a half-century. When a OB population 
can be studied spectroscopically and individual stellar absorptions are 
found, then main sequence fitting provides a reliable distance. But when 
they are too obscured or unavailable, then distance estimation methods 
are varied, unreliable and inconsistent. Consider, for example, the YSC 
surrounding the Herbig AeBe star LkHa 101, one of the brightest in- 
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Figure 2. The universal cluster X-ray luminosity function (XLF) found in Chandra 
YSC studies. NGC 1333 (77 stars from Getman et al. 2002), IC 348 (168 stars from 
Preibisch & Zinnecker 2002), and the ONC (1508 stars from COUP data truncated 
at \ogLx ~ 31.5 erg/s, Getman et al. 2004). 



frared sources in the sky, which illuminates the nebula NGC 1579. The 
cluster distance has variously been estimated to be 140 pc, >800 pc, 
~ 340 pc, and ~ 700 pc (Tuthill et al. 2002, Herbig et al. 2004). If it 
lies at the nearer end of this range, it is one of the closest YSCs. 

The X-ray luminosity function (XLF) of YSCs has two remarkable 
empirical characteristics that should render it an effective and accurate 
distance estimator for clusters such as this. First, the shapes of differ- 
ent YSC XLFs appear to be remarkable similar to each other, once a 
richness-linked tail of high luminosity O stars is omitted (Figure 2). Sec- 
ond, the shape of this 'universal' XLF strongly resembles a lognormal 
with mean < logL^; >~ 29.5 erg/s (0.5 — 8 keV band) and standard 
deviation ^(logL^;) ~ 0.9 (Figure 3). 

While we do not have astrophysical explanations for either of these 
YSC XLF properties, the shape of the XLF can be roughly understood as 
a convolution of the IMF which breaks from the Salpeter powerlaw below 
~ 0.5 Mq (and itself is sometimes modeled with a lognormal curve) and 
the correlation between and mass. These two effects result in a 
steep falloff in the number of fainter X-ray stars in a YSC. This was 
dramatically demonstrated in the Chandra studies of the ONC where a 
factor of 10 increase in limiting sensitivity from logL2,(lim) ~ 28.0 erg/s 
in the early observations to ~ 27.0 erg/s in the deep COUP observation 
led to only a very small rise in detected lightly absorbed ONC stars. 
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Figure 3. The Orion Nebula Cluster (ONC) XLF with lognormal fit given in the 
text. COUP sample of 1528 sources with absorption corrected 0.5-8 keV luminosities 
(Getman et al. 2004). All X-ray luminosities are in the 0.5 — 8 keV band corrected 
for absorption measured from the X-ray spectrum. 

The peak of the XLF at < log >^ 29.5 erg/s can be used as a stan- 
dard candle for distance estimation, much as the lognormal distribution 
of globular cluster optical luminosities is used as an extragalactic dis- 
tance estimator (Harris 1991). One counts cluster members in flux units 
(this is the X-ray log A^— log S curve), matches the observed log A^— log S 
to the universal XLF, and reads off the distance from the offset. The 
table below gives a simple example for a hypothetical cluster at various 
assumed distances observed with Chandra for 100 ks. The tabulated 
values are the fraction of sources seen in broad count rate bins assuming 
typical PMS X-ray spectra and negligible absorption. 



Table 1. Simulated X-ray source counts for young stellar clusters 



logia;(lim) (erg/s) 


28.0 


28.5 


29.0 


29.5 


30.0 


Distance (pc) 


400 


710 


1270 


2250 


4000 


Flux bin (cts) 






Percent of sources 






5-50 


22 


46 


60 


76 


86 


50 - 500 


46 


42 


34 


22 


14 


500 - 5000 


28 


12 


6 


2 





5000 - 50000 


4 















The table shows that the observed logA^ — logS shape will differ 
dramatically for YSCs at different distances. For close distances or un- 
usually long exposures (as in COUP), most cluster members are quite 
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bright and few faint ones are seen. For far distances or short exposures, 
most cluster members are near the detection hmit which is 3 — 10 pho- 
tons for typical Chandra fields. The principal challenges are corrections 
for absorption, which can be derived (except for the faintest sources) 
for individual sources from the X-ray spectrum, and the elimination of 
extragalactic background sources, which can be achieved by the absence 
of a stellar counterpart in sensitive i^-band images. Accurate distances 
and error analysis would not be performed in broad flux bins as in the 
table above, but would be based on unbinned sources fluxes using a max- 
imum likelihood method as described by Hanes k, Whittaker (1987) and 
Cohen (1991). 

Acknowledgements: We thank Thomas Preibisch (MPIfR) for use of 
unpublished COUP results. This work was supported by NASA contract 
NAS8-38252 (Garmire, PI) and COUP grant SAO GO3-4009A (Feigel- 
son, PI). 

References 

Cohen, A. C. 1991, Truncated and censored samples: theory and applications. New 

York: M. Dekker 
Feigelson, E. D. & Montmerle, T. 1999, ARAA, 37, 363 

Feigclson, E. 2003, in Stars as Suns: Activity, Evolution and Planets, lAU Symposium 

219, 27 

Feigelson, E. D., Gaifney, J. A., Garmire, G., Hillenbrand, L. A., & Townsley, L. 2003, 
ApJ, 584, 911 

Feigelson, E. D. & Lawson, W. A. 2004, ApJ, in press, (astro-ph/0406529) 
Flaccomio, E., Damiani, F., Micela, G., Sciortino, S., Harnden, F. R., Murray, S. S., 

& Wolk, S. J. 2003, ApJ, 582, 398 
Flaccomio, E., Micela, G., & Sciortino, S. 2003, A&A, 397, 611 
Getman, K. V., Feigelson, E. D., Townsley, L., Bally, J., Lada, C. J., & Reipurth, B. 

2002, ApJ, 575, 354 
Getman, K. V. and 23 coauthors, 2004, ApJS, accepted 
Hanes, D. A. & Whittaker, D. G. 1987, AJ, 94, 906 
Hanson, M. M., Howarth, I. D., & Conti, P. S. 1997, ApJ, 489, 698 
Harris, W. E. 1991, ARA&A, 29, 543 

Herbig, G. H., Andrews, S. M., & Dahm, S. E. 2004, AJ, 128, 1233 

Jiang, Z., et al. 2002, ApJ, 577, 245 

Lada, E. A. & Lada, C. J. 1995, ApJ, 109, 1682 

Slesnick, C. L., Hillenbrand, L. A., & Carpenter, J. M. 2004, ApJ, 610, 1045 
Townsley, L. K., Feigelson, E. D., Montmerle, T., Broos, P. S., Chu, Y., & Garmire, 

G. P. 2003, ApJ, 593, 874 
Tuthill, P. G., Monnier, J. D., Danchi, W. C., Hale, D. D. S., & Townes, C. H. 2002, 

ApJ, 577, 826 



